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^ ' Abstract: Resummation and power-corrections play a crucial role in the plie- 

rs I nomenology of event-shape variables like the thrust T. Previous investigations 

showed that the perturbative contribution to the average thrust is dominated by 
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O I gluons of small invariant mass, of the order of 10% of Q, where Q is the cen- 

1-^ I ter-of-mass energy. The effect of soft gluons is also important, leading to a 

I I non-perturbative 1/Q correction. These conclusions are based on renormalon anal- 

ysis in the single dressed gluon (SDG) approximation. Here we analyze higher mo- 
^ ' ments of the thrust distribution using a similar technique. We find that the char- 

acteristic gluon invariant mass contributing to <(1 — T)™> increases with m. Yet, 
^ ' for m = 2 this scale is quite low, around 27% of Q, and therefore renormalon resum- 

d I mation is still very important. On the other hand, the power-correction to <(1— T)^> 

from a single soft gluon emission is found to be highly suppressed: it scales as l/Q^- 
In practice, <(1 — T)^> and higher moments depend also on soft gluon emission from 
configurations of three hard partons, which may lead to as{Q'^)/Q power-corrections. 
This issue is yet to be investigated. 
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1. Introduction 

The systematic analysis of event-shape variables in e~^e~ annihilation has become an 
active research field in the recent years . High quality experimental data in a wide 
range of center-of-mass energies Q (from 10 up to 200 GeV) provide an opportunity 
for precision test of the theory. On the theoretical side, the simple set-up of e~^e~ 
annihilation allows to probe directly the QCD vacuum and learn about the interface 
between perturbative and non-perturbative physics. 

Empirically, it is known for quite some time that perturbative approximations 
to moments of event-shape variables are consistent with experimental data only upon 
including additive power-corrections. A classical example is the case of the thrust, 
defined by 

T = 5l^4^ (1.1) 

Ei \Pi\ 

where the summation is over all the final-state particles Pi. The thrust axis ut is a 
unit vector which is set, for a given event, such that T is maximized. Below we shall 
use the variable t = 1 —T which vanishes in the two-jet limit. 

For the average thrust, as for several other average event-shape variables, one 
finds 1/Q corrections, <t> ~ <t>pT -l- X/Q, where A is some hadronic scale. These 
corrections are quite large at LEP energies. The case of average event-shape variables 
was intensively studied both theoretically Pj-fl^ and experimentally. 

Here we investigate perturbative and non-perturbative aspects of higher moments 
of 1 — thrust, <t™>. This subject was addressed before (see e.g. I]!!], [1^, |13|, p!^), 
but it did not receive yet the appropriate attention. Whereas some experimental 
results are already available |T5| , p!6[] , many more can be extracted from stored LEP 
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data |T^. This analysis is worthwhile: alongside interesting theoretical insights, 
the investigation of higher moments of event-shape variables may lead to additional 
precision measurements of as- 

Existence of a well-defined perturbative approximation to event-shape variables 
is guaranteed since they are infrared and coUinear safe. However, the remaining sen- 
sitivity of these observables to soft and collinear emission is usually quite high and it 
shows up in the form of large perturbative coefficients. For event-shape distributions 
close to the two-jet limit (the Sudakov region) the leading logarithms in the perturba- 
tive coefficients can be resummed to all-orders ||18|. Out of this special part of phase- 
space, and in particular when moments of event-shape variables are considered, only 



the leading and next-to- leading perturbative terms are known |T9| , |20| , pT| . Since the 
apparent convergence of these perturbative expansions in a standard renormaliz at ion- 
scheme and scale is slow, resummation seems necessary. A relevant source of large 
coefficients are renormalon diagrams [^, |^ |T^, ^ which reflect the effect of 
the running-coupling. 

One can imagine |]^, |2^, in analogy with the skeleton expansion in the Abelian 
theory, a possibility of reorganizing the perturbative expansion such that all diagrams 
which correspond to dressing a single gluon are summed flrst, then diagrams with 
two dressed gluons, etc. Formally, a systematic expansion of this type has not yet 
been shown to exist. However, the spirit of the Brodsky-Lepage-Mackenzie (BLM) 
scale setting method one can attempt to identify and resum the single dressed 
gluon (SDG) terms to all orders. This is the basis of the renormalon resummation 
methodology, which was applied in many different QCD examples [0] . In the context 
of event-shape variables in e^e~ annihilation we mention the case of the longitudinal 
cross-section ||2^ and that of the average thrust 0. 

Whatever resummation procedure is applied, it is clear that perturbation theory 
alone cannot predict the observed values of event-shape variables: the perturbative 
calculation uses quark and gluon flelds while the measurement is of hadrons. The 
effect of hadronization on the observables can be modeled by Monte-Carlo simula- 
tions. However, in order to gain some understanding on the nature of conflnement in 
QCD it is favorable to analyze raw hadronic data, and parameterize non-perturbative 
effects in the simplest possible way. 

It was noticed that hadronization does not involve large momentum transfer. 
Therefore perturbative results calculated in terms of partons can almost be directly 
compared with the data 0. Due to the sensitivity of the observables to soft physics 
some modiflcation of the perturbative result is still necessary. At the perturbative 
level, infrared sensitivity shows up in the form of infrared renormalons fl^: the 
coefficients increase fast, asymptotically as n!, and have a constant sign pattern. 
Consequently, the series is non-summable. Different resummation prescriptions, or 
"regularizations" of the perturbative sum, differ by power-suppressed terms. This 
ambiguity must be compensated at the non-perturbative level. In the SDG ap- 
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proximation, one can calculate the perturbative sum and extract the form of its 
ambiguity. This way 0, ^ a perturbative calculation can be used to identify the 
parametric dependence of non-perturbative corrections on the external scale Q. The 
magnitude of these power-corrections cannot be calculated, but assuming some uni- 
versality properties |^ ^ it is possible to estimate it for a large class of observables 
based on experimental data for one of them. 

The SDG approximation, as applied in in the case of the average thrust, 
provides a systematic framework for the analysis of perturbative running-coupling 
effects together with the related power-corrections. These two aspects of improving 
the truncated perturbative result naturally complement each other: they reflect the 
same physical phenomenon. 

It is clear that any observable, including the moments of 1 — thrust considered 
here, may eventually depend on other non-perturbative effects, which are unrelated 
to ambiguities of perturbation theory. We assume that these effects are not impor- 
tant. Since the power-corrections identified from the ambiguities of the perturbative 
expansion have a definite Q dependence, this assumption can be confronted with 
experimental data. 

It should be emphasized that, contrary to the Operator Product Expansion 
(OPE), the SDG renormalon approach lacks the rigor of a systematic expansion: 
there is no small parameter which distinguishes the contribution of multiple emis- 
sion from that of a single emission. One should be aware of the possibility that in 
certain cases the leading power-corrections cannot be analyzed in the framework of 
the SDG calculation as they are associated with soft gluon emission from configura- 
tions of three of more hard partons. In particular, in the case of the higher moments 
of 1 — thrust, such effects may yield Ois{Q'^)/Q power-corrections |]14|, |, |Tl|, |T2|. We 
shall return to this important issue in the conclusions. 

Further subtleties are related to the fact that event-shape variables are not com- 
pletely inclusive 0: they are sensitive to certain details of the final state, while 
the resummation procedure we use is completely inclusive with respect to the frag- 
mentation products of the gluon. 

The purpose of this work is to examine, in the framework of the SDG approx- 
imation, running-coupling effects and power-corrections to higher moments of 1 — 
thrust. We begin, in section 2 by analyzing the thrust distribution. As in [0, where 
the average thrust was analyzed, we use the "massive" gluon dispersive approach 
[^]. In this framework resummation as well as parameterization of power-corrections 
are obtained from a so-called characteristic function. After explaining the main as- 
sumptions of the SDG model we calculate the characteristic function for the thrust 
distribution. We then devote a short discussion to the domain of applicability of the 
SDG result as a function of the thrust. Next, in section 3 we evaluate the charac- 
teristic functions for various moments of 1 — thrust and study their properties. In 
particular we extract the characteristic invariant mass of gluons contributing to the 
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various moments and compare the significance of running-coupling effects to other 
perturbative contributions. We also quantify the effect of the non-inclusive contribu- 
tion at the next-to-leading order. Finally, we extract the leading power-corrections 
implied for the various moments of 1 — thrust by the ambiguity of the perturbative 
SDG result, and identify the regions of phase-space from which they emerge. The 
conclusions are given in section 4. 



2. The thrust distribution in the single dressed gluon model 



The differential cross section with respect to the thrust da/dt can be calculated in 
perturbation theory in the single dressed gluon (SDG) approximation as follows 

/o e Jo e 



da , , 



Cf /' - K.W^)T((, t) =Cf t- p(cQ'') [^(t, t) - r(0, t)] (2.1) 

Jo e Jo e 



SDG 

2 /n2 



where the integration over e = /i /Q corresponds to inclusive summation over final 
states into which the emitted gluon fragments. In this calculation a final state 
of invariant mass fi"^ is characterized by a single universal (i.e. observable indepen- 
dent) function p(/i^) which is identified as the time-like discontinuity of the coupling 
a[k'^) = as[k'^)/7T, where the bar stands for a specific renormalization-scheme. The 
two integrals in (p.l|) are related by integration by parts: 

^(e,t) = -e^^(e,t), (2.2) 



and the "time-like coupling" a^gi^^"^) obeys 

/i , 2 = p(/i ). (2.3) 



The thrust characteristic function JF(e, t) is obtained from the following integral 
over phase-space, 

T{e,t)=[ dxidx2M{xuX2,e)6{l-T{xi,X2,e) -t) (2.4) 

v/phase space 

where CpaAi is the squared tree level matrix element for the production of quark- 
anti-quark pair and a gluon of virtuality /i^ = eQ^, and 



M{xi,X2,e) = ^ 



(xi + e)2 + (x2 + e; 
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[2.5) 



(1-Xi)(l-X2) (1-Xi)2 {1-X2)\ 

The integration variables Xi^2 represent the energy fraction of each of the quarks in 
the center-of-mass frame. The energy fraction of the gluon is X3 = 2 — xi — 0:2. 

For the calculation of the characteristic function we shall use the following 
definition of the thrust 



(2.6) 
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In case of three partons (a quark, an anti-quark and a "massive" gluon) it yields [0], 
1 - T(xi, X2, e) = min |l - xi , I - X2 , 1 - (2 - xi - X2y - 4e| . (2.7) 
Note that in the definition ( |2.(j| ) the denominator is modified with respect 



to the standard one (p..l|), Y.i — ^ Y^iEi, in a way that does not change the 
observable for massless partons. The virtual ("massive") gluon is understood to 
fragment eventually into massless partons. This modification ensures that the value 
of the thrust calculated with a "massive" gluon will be correct, provided that all 
the (massless) partons produced in the process of the gluon fragmentation end up 
in the same hemisphere with respect to ray The inclusive calculation performed 
here is justified only if the gluon fragmentation is predominantly coUinear, in which 
case fragmentation into opposite hemispheres will be rare. The discrepancy between 
the inclusive "massive gluon" calculation and the full non-inclusive calculation was 



addressed before by several authors |]^, |^ |T^ . This issue will be further discussed 
in the next section in the context of the moments of 1— thrust. 

The last ingredient for the calculation of the characteristic function is the phase- 
space^. Fig. shows the three parton phase-space in case of a gluon with a "mass" of 
fi"^ = eQ^ = 0.1 Q^. The external boundaries of phase-space in the figure correspond 
to the softest gluons (the upper curved line) and to the hardest ones (the lower 
linear line). The dashed lines represent the separation of phase-space according to 
which particle carries the largest momentum and thus determines the thrust axis (cf. 
eq. (p.Tp): in the upper left region T = X2, in the upper right region T = Xi and in 
the lower region T = \J x^ — 4e. 

It is important to note that the two upper regions of phase-space in the figure 
(where one of the primary quarks carries the largest momentum, T = Xi 2) have 
three corners with a definite physical meaning for e <^ 1: 

(i) the collinear limit corresponding to a two-jet configuration, with t = t. This is 
the lowest possible value of t, given e. 

(ii) the three parton symmetric limit, corresponding to a three-jet configuration, 
t = (2^/1 + 3e —1^/3. This is the highest possible value of t, given e. 

(iii) the large- angle soft gluon limit, with t = -^e. Here gluons are emitted close to 
90 degrees with respect to the quark-ant i- quark direction. 



It is clear from fig. |I] that in order to perform the integral in ( ^.41 ) one has to treat 
separately values of t smaller than y/e where the soft gluon phase-space boundary 
(the curved line in fig. |l|) is relevant, vs. t larger than ^/e, where it is not. As 
explained above, the value t ~ ^/e corresponds, for small e, to the situation where 

^The reader is referred to m for more details on the calculation of the phase-space boundaries. 
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Figure 1: Phase-space for the emission of a "massive" gluon with = eQ^ = 0.1 
in the plane of the quark and anti-quark energy fractions (a;i^2)- The outer continuous 
lines represent phase-space hmits, the dashed hues separate the phase-space according to 
the identification of the thrust axis nx and internal continuous lines are constant thrust 
contours. These contours are drawn with separation of AT = 0.03, stating at the lowest 
possible value ofl — T = e = 0.1 where two of the three partons are roughly collinear, up 
to the highest possible value of 1 — T = {2\Jl + 3e ~ /3 — 0.427 where the momentum 
splits equally between the three partons. 



the gluon is soft and emitted in a large angle with respect to the quark-ant i- quark 
direction. This part of phase-space is the source of the 1/Q power-corrections for the 
average thrust, and as we shall see in the next section, also for the leading, though 
much suppressed, infrared power-corrections for higher moments of 1 — thrust, at 
the SDG level. Note that other limits of phase-space, as well as the squared matrix 



element (2^) itself do not contain any y/e terms. In particular, in the collinear limit 
t ~ e. 

Evaluating (|2.4| ) we obtain the characteristic function for the thrust distribution, 
where the dominant contribution J^Q{e,t) corresponds to the phase-space regions 
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where one of the primary quarks carries the largest momentum (T = Xi^2) and 
J-'G{e,t) corresponds to the region where the gluon momentum is the largest (see 
fig. p. The superscripts / and h on J^q^e, t) denote low and high t values, respectively. 
These functions are given by 



[l-t + e? + (l + e? 



t ( gill 

In + 3-2- + -- + -t-g 

q-t \ t It 2 ^ 




t (g - 1) V t 2 t ^ 




l-t 



1 + e)2 + (1 + e - qf) In ^— ^ + (2t - g) ( g + 7 + 



(2.9) 



t 



t q-t 



where q = ^ {1 — ty + 4,e. 

Eq. ( |2.1| ) with jF(e, t) given by ( pl8| ) and ( |2.9| ) can now be used to calculate the 
SDG contribution to the thrust distribution. The leading order perturbative result 
p7[] can be recovered from (|2.1|) by assuming a constant coupling a^j((/i^) ~ const. 
This yields 



da , , 



LO 



/•I . 

CFacft(yU^) / — J^(e,t) 
JO e 



with 



^(0,t) 



- 3t 



2 ml 



2t 



C^a,ft(/i2)^(0,t) 



3(l + t)(l-3t) 
2 



(2.10) 



(2.11) 



t{l-t) t 2 t 

for t < 1/3. Eq. (|2.10| ) is the contribution to the thrust distribution of an on-shell 



gluon. The improved calculation ( |2.1| ) takes into account the emission of "massive" 
gluons that later dissociate. For a given thrust value, the gluon virtuality /i^ = eQ"^ 
can change in the following range: 



t < 
t > 



0<e<t 



f + it - i < e < t 



(2.12) 



Whereas the leading order result ( |2.1CI| ) depends on an a priori arbitrary choice of 
the renormalization scale /x^, eq. ( p^.l| ) resums running coupling effects to all orders, 
and is therefore renormalization-group invariant"'". Note that performing the integral 



■'"In fact, to have full renormalization-group invariance we have to fix the scheme of a. This 
requires further assumption concerning the diagrams that contribute to the gluon fragmentation 
("dressing the gluon" ^). Belo"w we fix it according to the Abelian-limit, in order to guarantee 
the correct resummation of the terms leading in Pq. 
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in eq. ( |2.1|) can be in turn re-formulated as a specific choice of the renormahzation 
point /i in ( p.lO| ) according to the BLM criterion ||25|. For the thrust distribution 



this representation seems rather cumbersome, since the BLM scale will depend on t 
in a complicated way. We shall use the BLM formulation later on for the moments 
of 1 — thrust. 

In spite of the resummation achieved by eq. ( p. ID , since this calculation takes 
into account just a single gluon emission its applicability is limited as follows: 

(i) At "high" values 1 — thrust, t ^ 1/3, only non-planar configurations, e.g. more 

than three hadron jets, contribute. The leading-order perturbative result ( p.lO| ) 
corresponding to three partons vanishes identically for t > 1/3 and the first 
non- vanishing contribution appears at the next order (9(a^). The resummation 
result (|2.1| ) does not vanish above t = 1/3. This fits the intuitive expectation: 
the fragmentation of the "massive" gluon is not restricted to the three parton 
plane. On the other hand, it is clear that in this region emission of two hard 
gluons from the primary quarks is important. 

(ii) At low values of 1 — thrust, Cpa lii^(l/^) ~ 1, the contribution of multiple 
emission of gluons which are both soft and coUinear becomes dominant and 
must be resummed to all-orders in perturbation theory [|18|. At even lower 



1 — thrust values, f3oa ln(l/t) ~ 1, this resummation breaks down since non- 
perturbative soft gluon emission becomes important. 

In view of these facts, we can expect the SDG approximation to describe the thrust 
distribution only in some restricted range of intermediate thrust values. Since soft 



and coUinear gluon resummation |T8| is complementary to dressing the gluon ( p.l| ), 
the two resummation procedures can be combined to describe the thrust distribu- 
tion in a wider range. The most promising approach to describe the differential 
distribution of event-shape variables in the Sudakov region (here small t) is based 
on introducing a non-perturbative "shape-function" (SF) |2^, |10|, |ll|, [1^. The phys- 
ical distribution is then obtained by convoluting this function with the perturbative 
result. Further analysis of the thrust distribution along these lines will appear in a 
separate publication. Here we proceed by analyzing moments of 1 — thrust. 



3. Moments of 1 — thrust 

Resummation of running-coupling effects in the SDG approximation was already 
demonstrated to be important for the average thrust. In particular, it was shown 
in that the resummation modifies significantly the value of as extracted from 
experimental data. The application of the SDG approximation to the average thrust 
was justified by asserting that low t = 1 — T values are suppressed in the average, 
and therefore the effect of multiple soft and coUinear gluons should not be significant. 
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From this respect the calculation of higher moments of 1 — thrust, <t™>, based on 
a SDG is even safer since they have a stronger suppression of the low t region. 
On the other hand, the contribution to <t™> from extremely high values of 1 — 
thrust, if: ^ 1/3, becomes more significant as m increases. For high enough m this 
contribution must be dominant, making the SDG calculation unreliable. 

Let us now assume that the contribution to <t"^> is dominated by intermediate 
values of 1 — thrust where the SDG approximation to the distribution (eq. ( |2.1D ) 
holds, and use it to calculate the m-th moment of 1 — thrust according to 



da 
'•1 de 



<r> ~ <t"^>sDG = f^{t) t'" dt = CF [ r dt C - a.«(eQ2)^(e, t) 

J dt J Jo e 



Cf / -aoff(eQ')-^<t™>(e), (3.1) 
JO e 



where in the last step we changed the order of integration defining 

T<t^> = J J^{e,t)t"'dt. (3.2) 



In order to use the resummation formula ( p.l[ ) we have to specify the time-like 
coupling dag{n'^). As explained in ||^, the renormalization-scheme defining the 
coupling d should be uniquely determined once an Abelian-like skeleton expansion 
is established in QCD. For the purpose of the current investigation we shall just use 
the one-loop approximation of the space-like coupling 

Po In [k^/A^j 

with A set such that d coincides with the Gell-Mann Low effective charge in the 
Abelian (large /?o) limit. This is realized, for example, if d is related to the MS 
coupling by a scale-shift, i.e. 

The corresponding time-like coupling is given by 



Po 



2 ' 



2 TT U ^A2 



(3.5) 



Using this coupling in ( p.l|) , the terms which are leading in Pq will be resummed 
correctly (ignoring the non- inclusive nature of the observable 0, which will be 
discussed below), while other terms which are sub-leading in jSo will be neglected. 

Before proceeding with the SDG analysis it is worthwhile to examine the known 
next-to-leading order result [|19|. The next-to-leading order coefficient is calculated 
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based on numerical integration over the three and four parton phase-space 
In MS it is given by 0, || 



<t> = 0.7888 Cf Ams (Q^) + (-1.1570 C/ + 4.4708 Cp Ca - 0.8445 Cp Nj') a|s (Q^) 
<t2> = 0.0713 Oms (Q^) + (0.3073 + 0.3280 Cp Ca - 0.0583 Cp A^/) (Q^) 
<t3> = 0.0112(:7^aMs (Q^) + (0.06877 + 0.04973 C^C^- 0.00808 C^A^/) a^g (Q^ 



<t4> = 0.0022 CpaMs((5 )+ 0.01622 + 0.00989 C^C^- 0.00145 C^A^/ a^siQ 



(3.6) 



In order to correctly identify the terms that originate from the running- coupling ^ 
we use ( p.4| ) to translate ( p.6|) to the a scheme. Then, expressing the Nf dependence 
of the next-to-leading order coefficients in terms of Pq = (11/12)C^ — (1/6)A^/ as 
done in the "naive non-Abelianization" procedure we obtain 



<t> 



<t'> 



Cp 
Cp 
Cp 
Cp 



0.7888 a(g2) + (3.7526 (3q - 1.1567 Cp - 0.1740 Ca) a^{Q^] 
0.0713 a(g2) + (0.2308 Po + 0.3073 Cp + 0.00762 Ca) a^iQ^] 
0.0112 a(g2) + (0.02981 Pq + 0.06877 + 0.005271 Ca) c?{Q'^\ 
0.0022 a(Q2) + (0.005045 /5o + 0.01622 Cp + 0.00191 Ca) a^iQ^, 



(3.7) 



Let us now look at the relative magnitudes of the different terms in the next-to- 
leading order coefficients in (|3.7|). Substituting the QCD values of Cp, Ca and Pq 
(for Nf = 5) we obtain the numerical values summarized in table ||. We find that the 





<t> 


<t2> 






Po 


7.1925 


0.4423 


0.0571 


0.00967 


Cp 


-1.542 


0.4098 


0.0917 


0.02163 


Ca 


-0.522 


0.0229 


0.0158 


0.00574 



Table 1: Decomposition of the next-to-leading order coefficient multiplying Cpa'^{Q'^) in 
eq. ( |3.7| ) for each of the moments of t = 1 — T. Here we substituted Ca = 3, Cp = 4/3 
and Nf = 5. 

non-Abelian CpCa component is always small. The relative significance of the CpPq 
term associated with the running- coupling decreases with m: this term is absolutely 
dominant in the case of the average thrust. It is still the largest in the case of <t^>, 
but it is no longer dominant. This trend continues for higher moments <t™'>, m = 3 
and 4. On the other hand, the significance of the double emission term, proportional 
to Cp^, increases with m. This fits the intuitive expectation: higher moments of 1 — 
thrust become sensitive to spherical configurations, and in particular, to multi-jet 
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configurations. This property, which makes the next-to-leading order perturbative 
series less rehable as m increases, also implies that the significance of the resummation 
of running-coupling effects by the SDG formula ( |3.1| ) compared to other contributions 
decreases with m. Note, however, that with increasing orders in perturbation theory, 
the terms associated with the running-coupling are expected in general to increase 
fast and eventually dominate the coefficients. 

We now want to calculate the characteristic functions (|3.2|) for the first few 
moments. Using eq. (|2.8|) we have 



<t"^>[^) 



(3.8) 



or, after rearranging the terms. 



where the difference 



J-G(e,t)+^^(e,t) 



r dt- 



t"" dt 
(3.9) 



•^Q(e,t)--^Q(e,t) = - 



t 



'l-t + e? + (l + e? 



ln-+2+-+e 2- 



t2 



e 



-^2 + 



2t 



contains the essential information on the region t ~ ^/e. As we shall see below, 
non-analytic terms in the expansion of J^<t'">(e) at small e imply power-suppressed 
contributions to <t™>. It is the ^/e in the upper boundary in the second term 
of (|3.9| ), corresponding to large-angle soft-gluon emission, which is the source the 
leading non-analytic terms. 

The resulting characteristic functions for m = 1 and 2 are shown in fig. ^ and 
these for m = 3 and 4 - in fig. 0. On this basis, the SDG perturbative sum in ( |3.1| ) 
can be calculated. As an alternative to evaluating the integral we can expand the 
coupling d^ff{eQ^) under the integration sign, e.g. in terms for a (Q^), obtaining. 



do a{Q') + di/5o <x\Q'') + ( d2 - :-do ) (31 a\Q'') + ■ ■ 



(3.10) 



where the coefficients are expressed in terms of the log-moments of the corresponding 
characteristic function, c/q = ^<t'">(0) and 



di 



:-lne)^ ^. 



(3.11) 



The log-moments can be evaluated by a straightforward numerical integration. The 
resulting values are summarized in table |[ In the cases of the average thrust and 
<t^> one identifies already at the first few orders the characteristics of infrared 
renormalons, namely fast increase (that asymptotically becomes factorial) and a 
constant sign pattern. This behavior is much delayed for the higher moments. 
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Figure 2: The derivative of the characteristic functions J-{e) for <1 — T> (upper plot) and 
<{l-Tf> (lower plot) function of logxo(e), where = eQ^ is the gluon virtuality. 
The dashed lines in the two plots represent the 0(e^/^) and ©(e^/^) approximations to 
^(e), respectively. 





<t> 








do 


0.7888 


0.0713 


0.0112 


0.0022 


di 


3.588 


0.1875 


0.0160 


0.000139 


d2 


21.06 


0.6296 


0.0218 


-0.00903 


ds 


154.03 


2.5922 


-0.0020 


-0.0555 


di 


1368.38 


12.721 


-0.3254 


-0.3108 


d5 


14464.59 


72.835 


-3.1238 


-1.8546 



Table 2: First few log- moment of the characteristic functions for various moments of 
t = l-T. 



The approximation of <t"^> by <t"*>sDG in ( |3.1| ), or by (|3.10|) , can be im- 
proved in a straightforward manner by matching it with the next-to-leading order 
result ( p.7|) . Such a procedure was used in for the phenomenological analysis of 
the average thrust. The additional terms in ( p.7| ), correspond primarily to double 
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Figure 3: The derivative of the characteristic functions j-{e) for <(1 — T)^> (upper plot) 
and <(1— r)^> (lower plot) as a function of log]^Q(e), where ji^ = eQ^ is the gluon virtuality. 



gluon emission (a Cp^ term). They also include a non-Abelian CaCf component, 
which depends 0, on the identification of the coupling a beyond the Abelian 
limit, and finally, also some residual /3o dependent piece, which is related [0 to the 
non- inclusive nature of the thrust. 

Had the moments of the thrust been completely inclusive with respect to the 
fragmentation products of the gluon, the terms proportional to (3q should have been 
fully contained in the SDG result ( ^.l] ) or ( |3.1CI| ). As explained in section 5.2 in [0], 
in the inclusive case the log-moments di, for any i, are equal to the terms leading in 
(3q in the exact calculation. In the present case the next-to-leading order coefficient 
di calculated as a log-moment of the inclusive characteristic function differs from 
the actual /3o dependent term in (|3.7| ). The "massive gluon" inclusive treatment 
is justified only if the discrepancy between them is small. In the case of <t> this 
discrepancy was found [0] to be tiny. Comparing table | with eq. (|3.71) we find that 
it is only 4.4%! Evidently it increases for the higher moments of 1 — thrust. For 
<t^> the inclusive approximation still seems reasonable: the discrepancy is 18.7%. 
It is much worse for <t^>, where the discrepancy is 46% and it completely breaks 
down for <t^>. 
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The physical reason why the inclusive calculation for the high moments fails is 
that <t'"> becomes sensitive to large t values. The high t region is correlated with 
high gluon virtuality (see e.g. ( p.l2|) ). High virtuality allows the gluon fragmenta- 
tion products to spread in large angles compared to the original gluon momentum. 
Consequently, fragmentation into opposite hemispheres becomes more probable. In 
such full non-inclusive calculation would yield a higher thrust (T) value 

compared to the inclusive calculation that takes into account the gluon momentum 
itself (this follows from (|2.6|) using the triangle inequality). As a result, the inclusive 
calculation under-estimates the thrust distribution da/dt at large t on the expense 
of over-estimating it for somewhat lower t values. This is why for all the moments 
the full Po dependent term in ( p.7|) is larger than the inclusive result di. The success 
of the inclusive approximation for the average thrust and for <t^> is related to the 
fact that most of the gluon fragmentation is roughly collinear. 

For the first two moments of 1 — thrust, the coefficients in ( ^.10 ) increase fast 



due to infrared renormalons already at the first few orders. It is clear that a(Q^) 
is not a good expansion parameter. A simple way to approximate the perturba- 
tive sum (ignoring the renormalon ambiguity) is the BLM method. In this case we 
approximate ( p.lD by 

<t"'>sBG-CFd,s{fil^J —J^<t^>{e) = CFd,s{fjl^JJ^<tr^>{0) (3.12) 

JO e 

where, at leading order, the BLM scale is the center of the characteristic function, 
i.e. the average virtuality of the gluon contributing to <t"'>sDG) 

/ /•! . de / ■ de\ 

AtBLM = exp f - J^<jm>(e) (-Ine) — / -^<t™>(e)— 1 = Q'^ exp (-di/do) . 

(3.13) 

In principle, higher order log-moments ( p.ll|) can be used to improve the approxima- 
tion in ( p.l2| ) by further modifying the scale (see, for example, eq. (33) in p4|). For 
the qualitative discussion here, the leading-order BLM approximation will suffice. 
Using ( |3.4|) one can translate the BLM result to the MS scheme, 

<t™>SDG ^ C'fAmS (AtBLM,Ms)-^<t™>(0) (3.14) 

with — = ii'^ 



BLM, MS 



The first two log-moments in table ^ allow to calculate the leading-order BLM 
scales for the various moments of 1 — thrust according to ( p.l3|) . The results are 



shown in table |^. The numbers in the first row of the table have the interpretation 
of the average gluon virtuality contributing to <t"^>. We see that 

(i) the typical virtuality, which is also the correct argument for the coupling in 
( |3.12| ), is much lower than Q for the first moments of 1 — thrust, especially for 



the average thrust and for <t >. 
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<t> 


<t^> 


<t^> 


<t^> 


A'-blm/ Q 

/^BLM,MS / Q 


0.1028 
0.0447 


0.2685 
0.1167 


0.4889 
0.2125 


0.9689 
0.4211 



Table 3: BLM scales in the "skeleton scheme" and in MS for various moments oft = 1 — T. 



(ii) with increasing m highly virtual gluons become dominant (at large m, /Xblm — ^ Q)- 

These features can be read directly from the characteristic function curves. In fig. ^ 
we see that J^<t>(e) is wide and centered at very low virtualities while jF<i2>(e) is 
narrower and centered at higher virtualities. Fig. ^ shows that this trend persists 
for the higher moments. Note that J^<t^>{t) for m > 3, as opposed to m = 1 and 
2, is not positive definite. Thus for m > 3 the interpretation of /Xblm as the average 
gluon virtuality contributing to <V^> is no longer accurate. Note that for m = 4 
the cancelation between positive and negative contributions to di is already large, 
making the approximation of the integral ( |3.1| ) by the coupling at the BLM scale 
absolutely unreliable. 

Let us now turn to discuss the power-correction implied by the SDG perturbative 
sum. The starting point is the observation that eq. ( |3.1| ), understood as an all-order 
perturbative sum, is ill-defined due to infrared renormalons. The standard way 
to resum asymptotic series of this type is the Borel method. To obtain a Borel 
representation one expresses the coupling in (|3J]) as 



Klil^') = ^ r rf.a(.)^e-^^°^, (3.15) 
where the sin factor arises from the analytic continuation to the time-like region. For 



the one-loop coupling (|3.3|) , a{z) = 1. Substituting this in (|3.1|) and changing the 
order of integration one arrives at the following Borel-sum, 

<t"'>suG = -f dzB^tm^{z)e-'''-^7^ (3.16) 



/5o Jo 
with 

^ , , sin TTZ , • , . de , _s 

5<,™>(z) = / e-^<,™>(e)_. (3.17) 

TTz Jo e 

The integral over z in (|3.16|) is ill-defined yielding an ambiguous result (see 0, ^ 
for more details). The physical reason for this ambiguity is that our "perturba- 
tive" calculation ( p.l|) actually depends on the coupling at all scales, including the 
infrared. This ambiguity can thus be resolved only at the non-perturbative level. 
Nevertheless, in practice it is possible to cure the ambiguity by defining the per- 
turbative sum through some regularization procedure, like the principal-value or a 
momentum cutoff. At the next step, one adds to the regularized perturbative sum 
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explicit power-suppressed terms that have the same dependence on as the leading 
ambiguities and a normalization which is controlled by free parameters to be fixed 
by a fit. An attractive possibility is to write this parameterization in terms of 
the small moments of coupling d{k'^), assumed to be regular in the infrared at the 
non-perturbative level. Thanks to the assumed universality of a(fc^), this parameter- 
ization immediately implies universality of the power suppressed terms for different 
observables, up to calculable factors which depend on the characteristic function. 

Notice that since the time-like coupling ( p. 51) is finite for any /i^ and has an 
infrared fixed-point: acff(O) = l//3o, the integral (|3.1|) with (|3.5|) is well-defined. This 
is the so-called "Analytic Perturbation Theory" (APT) result. It differs by 



ambiguous power terms^ from the corresponding Borel-sum ( p.l6 ). The APT result 



just like the principal-value Borel-sum, is a specific regularization of the perturbative 
sum. Its advantage is that it is straightforward to calculate. 

In general, if we ignore sub-leading perturbative terms that are not included in 
the SDG approximation (|0|), the physical result should be given by the regularized 
perturbative sum plus power-corrections. The apparent ambiguity in the choice of 
the regularization procedure is eliminated by the power terms. In absence of relevant 
non-perturbative calculation, the only way to appreciate the significance of the power 
terms is by comparing different regularization procedures. 

One has to distinguish between two types of power-suppressed ambiguities of the 
perturbative sum: 

(i) differences between various regularizations, e.g. between the APT integral and 
the principal-value Borel-sum. 

(ii) the specific differences which are associated with infrared scales, i.e. space-like 
momentum scales at which the coupling is not under control in perturbation 
theory. 

This second type of ambiguity can be studied by introducing a momentum cutoff 
in the space-like region 0, |^. It was shown to be related to non-analytic terms 
0, ^ H, 0, in the small gluon virtuality (or e) expansion of the characteristic 



function. In case of observables^ that admit an operator product expansion (OPE), 
the parametric dependence of these ambiguous terms on Q should fit the dimen- 
sions of the higher-twist operators. On the other hand, the first type of ambiguity 
may appear due to both analytic and non-analytic terms in the small e expansion 
of the characteristic function This implies that certain regularizations differ 



from others, as well as from the physical non-perturbative result, by non-infrared 
non-OPE power terms. The simplest possibility is that non-perturbative effects are 



§See 0, |g and refs. therein. 

^This does not apply to moments of 1 — thrust which cannot be expressed in terms of local 
operators. 
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restricted to large distances. Then a momentum cutoff in the space-hlce region as 
introduced in |3y, ^, or equivalently, the principal- value Borel-sum, define a class of 



favorable regularizations which differ from each other, as well as from the physical 
result, just by infrared power-corrections. We stress that the APT integral does not 
belong to this class of regularizations, and is suggestive of a different scenario for the 
regularization of perturbation theory, advocated in |^ . 



In order to study the ambiguities of perturbation theory that emerge from the 
SDG calculation in the case of moments of 1 — thrust (|3.1| ), we calculated the small e 
asymptotic expansions of the characteristic functions J^<t™>(e) given by eq. (|3.9| ). 
These expansions for m = 1 through 4 are summarized in the appendix in eq. (|A.lj) 



through ([A.4|) and in table ^. According to the explanation above, non-analytic terms 



in the expansion of the characteristic functions are of special interest, since they imply 
ambiguity of perturbation theory which is associated with large distances. There are 
two physically distinct sources of non-analytic terms contributing to the expansions 
( [A. 1|) through ( |A.4| ). The first is the contribution of coUinear soft gluon emission. 



It originates in the lower integration limit {t = e) of the two terms in ( ^.9| ). The 
second is large-angle soft-gluon emission, which originates in the upper integration 
limit (t = ^/e) of the second term in ( |3.9| ). To distinguish between these two regions 
of phase-space, we evaluated separately the second integral in (|3l^ ) at its upper limit. 
The results, which represent the contribution of large-angle soft-gluon emission to 
^<t™>(e), are summarized in eq. (|A.5|) . The expansions in the appendix can be used. 



based on the general formulae obtained in ref. [^, to analyze the power-corrections 
implied by 

Note first that all the leading non-analytic terms in the small e expansion of 
the characteristic functions are of the half-integer type. These terms, as well as all 
the other half-integer terms, originate exclusively in the large-angle soft-gluon region 
through the upper limit t = ^/e in the second integral in ( p.9| ). This conclusion 
follows from the comparison of eq. ( [A.5[ ) with the full expansion in eq. ( [A.l| ) through 



Note also that in eq. ( |A.1[ ) through ( [A.4| ) there are no logarithmic terms at 
order e In e. In general, presence of such terms in the characteristic function implies 
[^, ^ infrared power-corrections that scale as Their absence at the level of 

the moments of 1 — thrust, <t"^>, is not expected a priori, since they are present 
at the level of the thrust distribution, see eq. (^]8|) with ( p.9|) . Comparing eq. (|A.5| ) 
to eq. (^) it becomes apparent that for the average thrust the elne terms cancel 
between coUinear and large-angle soft terms. For higher moments of 1 — thrust such 
terms do not appear even in the separate contributions of these two phase-space 
regions (see eq. ( |A.5| )). The first non- vanishing logarithmic terms for the first two 
moments, m = 1 and 2, appear at order e^lne, namely l/Q^. These terms originate 
in both collinear and large-angle soft-gluon regions. The leading logarithmic terms 
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for higher moments, m = 3 and 4, appear much later, and exclusively" due to the 
collinear limit. Since in non of the moments, do logarithmic terms appear as the 
leading non-analytic terms we can safely ignore them in the more detailed analysis 
that follows. 

Let us now summarize some formulae [0 for the regularization dependence of 
using the simplest one- loop ansatz ( p.3|) for the coupling. Given a generic term 
of the form** Cne^ in the small e expansion of JF<tm>(e) — J^<t'">(0), the difference 
between the APT regularization and the Borel-sum is given by [^, |3y, 0] , 

5<r> = -c„-^ f^j exp(±mn). (3.18) 

Note that for half-integer values of n, this formula yields an ambiguous imaginary 
result. This is the ambiguity of the Borel-sum. The principal- value regularization 
amounts to taking the real part of the Borel sum, and thus it coincides with the APT 
integral, as far as these terms are concerned. On the other hand for integer values 
of n the difference between the APT integral and the Borel-sum is unambiguous and 
real. 

To quantify the ambiguities associated with infrared scales we define a space-like 
momentum cutoff with tj = In ^/ij/A^^. The contribution to the principal- value 
Borel-sum from momentum scales below fij is given by [0, 

^ Ci? / u? \ " sin vrn , , ^. , , ,^ ^ ^^ 

<t^>,^ = _c„ ^ exp{-nti) Ei{ntj). 3.19 

In accordance with the general statement above, <t'">iR is non-zero for non-analytic 
terms with half-integer n while it vanishes for analytic terms where n is integer. 

Using ( |3.18| ) and ( |3.19| ) with the expansions of the characteristic functions in 
eq. ( |A.1| ) through ( [A.4| ) we find the leading ambiguities associated with the pertur- 
bative sum ( ^.Ij ) of each moment of 1 — thrust. The two rows in table ^ correspond 
to ( p. 181 ) and ( |3.19| ), respectively. The table shows the parametric dependence of the 



leading ambiguity on the center-of-mass energy Q. Table ^ shows the corresponding 
numerical values of the power terms for Q = Mz, normalized by the (approximate) 
perturbative result ( p.l2| ) with the BLM scales of table |^, for yU/ = 2 GeV, assuming 
a,(Mz) =0.115 and Nj = 5. 

As mentioned above, in all the cases the leading contributions from the infrared 
( |3.19|) appear due to half-integer terms associated with large-angle soft-gluon emis- 



sion. These ambiguities imply existence of non-perturbative power terms which scale 
as l/Q for <t>, as for <t^> and <t^>, and as for <t'^>. 



"in eq. ( A. 5 ) for m = 3 and 4 there are no logarithmic terms at aU. 



"As stressed above we ignore logarithmic terms. Regularization of integrals like (3.1) in presence 
of such terms was also analyzed in [Q. 
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<t> 








(5<t™> 

<t™>IR 


i/Q 
i/Q 






i/Q' 



Table 4: Leading ambiguities in the perturbative calculation for various moments of 





<t> 


<t2> 






5<t™>/<t™> 

<t™>iR/<t™> 


±0.075 z 
0.17 


-0.3810-3 
0.1410-3 


0.2210-3 
-0.2410-3 


0.3710-3 
-0.35 10-6 



Table 5: Relative magnitudes of the ambiguities in the perturbative calculation for various 
moments oft = l — TatQ = Mz. 



Considering the leading ambiguity of the perturbative sum ( |3.18|) , a definite 



difference exists between the case of the average thrust, m = 1, and higher moments, 
m > 2. In the case of <t> the leading ambiguity of the Borel-sum is due to the 
same ^/e term which dominates the infrared contribution, whereas for <t'^> with 
m > 2 the leading ambiguity originates in an analytic term (cie) in J^<tm>(e). It is 
therefore not associated with large distance scales. Moreover, a careful examination 
of the source of the terms proportional to e, (compare e.g. ( |A.5| ) with the full 
expansion ( |A.1|) through (|A.4| )) shows that they are not associated with any definite 
part of phase-space. 

The numerical values in table ^ clearly indicate that all the power-corrections 
that appear in the SDG approximation are small for any m > 2. The relative error 
at Mz is less than a pro-mil. The case of the average thrust is unique having a 
significant contribution from the infrared: with the parameters quoted above it is 
17% (!) at Mz. We stress that in contrast with the average thrust case, the difference 
between the APT and Borel-sum regularizations for m > 2 becomes parametrically 
larger, as well as numerically larger at Mz, than the infrared contribution to the 
Borel-sum. This might give an opportunity to use experimental data in order to 
constrain l/Q"^ power terms such as the ones making the APT integral different 
from the principal-value Borel-sum. 



4. Conclusions 



We started in section 2 by calculating the characteristic function of the thrust dis- 
tribution, jF(e, t). The result, summarized by eq. 



rsD and (|2.9|) , forms the basis for 
resummation of running-coupling effects in the SDG approximation in this case. On 
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its own the SDG approximation is expected to describe the physical distribution only 
in a limited range of thrust values. However, when combined with other techniques 



IS, 12], the range of applicability of the calculated distribution can be extended. 



leading to more meaningful comparison with experimental data than available to- 
day. 

In section 3 we concentrated on the first few moments of 1 — thrust, <t™>, 
m = 1 through 4. The corresponding characteristic functions ^<t'">(e) were obtained 
by a straightforward integration of T[e^ t) and then used to study the properties of 
the SDG approximation to <t''">. We saw that the characteristic mass scale of 
gluons contributing to <t'"> increases fast with m (see table |]). Whereas typical 
gluon virtuality contributing to the average thrust is about 10% of the center-of-mass 
energy Q, this fraction becomes 27% for <t^>. Still, if one is using the MS scheme, 
the natural renormalization scale /Xblm.ms is quite far from the naive choice n = Q: 
for <t^> we have /Ublm.ms — 0.12(5- This means that the effect of the resummation 
compared to the naive perturbative treatment is still very significant in the case of 
<t'^>. 

For high moments m > 3, a significant source of uncertainty in the available 
next-to-leading order perturbative approximation (see eq. (|3.71 )) is related to multi- 
jet configurations that contribute at high values of t, t ^ 1/3. Improving the ap- 
proximation requires a full next-to-next-to-leading order calculation. 

Our renormalon analysis was performed in the inclusive "massive gluon" ap- 
proach. The advantage of this approach, apart from its simplicity, is that it naturally 
generalizes to include higher orders [0 in the f3 function of the "skeleton scheme". 
In addition, it allows for non-perturbative effects to be parametrized in a trans- 
parent way using an infrared regular running-coupling. Using this parameterization 
the universality assumption can be tested. Of course, the inclusive resummation is 
justified only if the non- inclusive effect is small. An alternative approach to perform 
renormalon resummation is based on the "naive non-Abelianization" procedure. One 
can define non-inclusive large Nf characteristic functions (see [0, for <t™>, and 
then restore the full /5o of the non-Abelian theory. Comparing the expansion of the 
inclusive resummation formula ( J!T) to the next-to-leading order perturbative expan- 
sion, we obtained a quantitative estimate of the non- inclusive effect. The discrepancy 
at the next-to-leading order is very small (4.4%) for <t> and it increases for higher 
moments. In the case of <t^> the inclusive approximation still seems reasonable (a 
discrepancy of 18.7%) but less so for <t^>. Due to the appreciable discrepancy it is 
worthwhile to compute also the non-inclusive resummation. 

We find that within the framework of the SDG calculation, power-corrections 
to <t™>, with m > 2, are highly suppressed (see tables ^ and The main ques- 
tion that remains open concerning the phenomenology of <t^> is the significance of 
power-corrections from configurations of three hard partons plus a soft gluon. Such 
corrections can be [^, [11], |T2|, |I^ as large as as{Q^)/Q- We further address this 
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subject below. 

In the case of the average thrust the leading ambiguity of the perturbative sum 
scales as l/Q. This ambiguity originates in a particular part of phase- space where a 
soft gluon is emitted at a large angle. Being associated with large distance physics, it 
is quite clear that this ambiguity can be resolved only at the non-perturbative level, 
by including explicit non-perturbative terms that fall as l/Q. We showed that for 
higher moments of 1 — thrust the infrared contribution to the SDG perturbative sum 
is dominated by the same large-angle soft-gluon region of phase-space. The corre- 
sponding non-analytic terms in the characteristic functions lead to much suppressed 
power-corrections, which scale as for <t^> and <t'^> and as for <t^>. 

For the various moments of 1 — thrust, <t™>, as for any other time-like ob- 
servable, there are further ambiguities in the summation of perturbation theory that 
are not related to large distance scales. Such ambiguities become apparent when 
comparing different regularizations, such as the principal-value Borel-sum and the 
APT integral. In the case of the moments of 1 — thrust, these ambiguities scale as 
l/Q^ and thus dominate over the infrared ambiguity for m > 2. We saw that this 
type of ambiguity is not associated to any particular part of phase-space. This sug- 
gest that, contrary to the soft emission ambiguities, it does not signal any genuine 
non-perturbative effects. It may be then possible to eliminate the ambiguity just 
by choosing the correct class of regularizations of the perturbative sum. The sim- 



plest possibility is that this class is defined by a space-like momentum cutoff [30 



or equivalently by the principal-value Borel-sum. One should be aware that other 
scenarios |^T| are possible as well. 

Finally, we would like to view our conclusions concerning the power-corrections 
for the moments of 1 — thrust in the SDG model in the context of previous analysis 
of the thrust distribution. For this purpose we briefly recall the results of refs. 
|ll], |12|. Assuming a two-jet configuration, it was shown that the main effect of non- 
perturbative soft gluon emission is a shift of the Sudakov-resummed perturbative 
spectrum to higher values of t\ 



da , ^ da 
dt^ ' dt 



(t-At), (4.1) 



where At = X/Q. A priori, this formula applies only in the range At <C t <^ 1/3. 
It was demonstrated [Q that the measured thrust distribution can be fitted over a 
large range of t and Q values by introducing such a shift. 

It is clear that ( [4. 1|) cannot describe the physical distribution at extremely small 
values of t, t ^ At, where multiple non-perturbative soft gluon emission becomes 
essential. This difficulty can be resolved ||Tl], |12[, by introducing a non-perturbative 
(observable dependent) shape-function (SF) to describe the energy flow in the final 
state. As explained in [0, in this case the physical distribution is obtained by a con- 
volution of the Sudakov-resummed perturbative spectrum with the shape-function. 
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The resulting distribution at extremely small t (to the left of the distribution peak) 
depends on the form of the shape-function but at higher t it approximately coincides 
with the shifted distribution ( [4.1|) . 

Both the simple shift model and the shape-function approach |Tl|, |T^, like the 



Sudakov-resummed perturbative spectrum on which they are based, strongly rely on 
the two-jet kinematics. These approaches strictly do not apply to the large t region 
where one gluon becomes hard. The success of the fits [Tl|, |12| in a large range of 
t and Q values is encouraging and it might suggest that also the first few moments 
of 1 — thrust could be studied in this framework. Yet, one should be aware of the 
fact that the distribution is rather fiat at large t (where one gluon becomes hard) 
and thus the total effect of the shift, or the convolution with the shape-function, is 
minor there. On the other hand the first moments crucially depend on this very same 
region of t. Therefore it is dangerous to analyze the non-perturbative corrections to 
the moments relying on the apparent success of the fits to the distribution. 

Nevertheless, it is interesting to see the consequences of the assumption that the 
simple shift model |j5| or the shape-function approach |jll], |12| apply beyond the small 
t region. In this case one finds |jTT|, |12|, [13| , 



<t>SF = <t>PT + \i/Q (4.2) 

<tSsF = <tSpT + 2Xi<t>p^/Q + A2/(5^ 

where the scales Aj are the moments of the shape-function and in the case of a 
shift ( ^4.1|) , they are simply Aj = A*. This can be compared to the infrared power- 
corrections found in the SDG calculation (table namely 

<t>suG = <t>PT + A/Q + 0{l/Q^) (4.3) 

<t'>SDG = <t'>PT + 0(l/Q'), 

baring the fact that the perturbative part <t™>pT can be quite different in the two 
approaches: in ( [4.2|) Sudakov resummation is implicitly assumed while in ([4.3|) SDG 
renormalon resummation is assumed. 

For the average thrust the two approaches predict the same type of leading 
power-correction. Let us examine the case of <t^>: 



the leading non-perturbative term in (|4.2|), 2Ai<t>pT/<5, is proportional to 



as{Q^)- Therefore, this term is attributed to soft emission around a config- 
uration of three hard partons. Clearly, this is beyond the scope the SDG 
calculation performed here. On the other hand, this term crucially depends 
on the extension of the distribution models of [^, [Tl|, |12| beyond their range 



of validity (two-jet kinematics). Therefore further theoretical work is required. 
Note that this prediction for the leading power-correction of <t^> is rather 
easy to verify, or exclude, using experimental data: due to the large ratio 
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(see eq. ( p^) ) between the normalization of <t>pT and <t^>pT in ([4.2|) , the 
numerical significance of this power-correction is quite large. 



(ii) the second term in ( [4.2|) , \2/Q'^, is not suppressed by any power of a^, and is 
therefore attributed to soft emission around the two-jet configuration. On the 
other hand, the SDG calculation ( ^I3| ) does not yield any infrared correc- 
tion to <t^>pT. There are two possibilities how these facts can be reconciled. 
The first is that once a more complete^^ description of the thrust distribution 
is achieved there will not be any infrared correction to <t^>. The sec- 

ond is that an infrared l/Q"^ correction will emerge entirely from double gluon 
emission. It is a possible, but yet an unusual, situation that double soft gluon 
emission becomes parametrically less suppressed than a single soft gluon emis- 
sion around the same hard configuration. Therefore it will be interesting to 
investigate power-corrections to <t^> from double gluon emission within the 
renormalon approach. 

We comment that techniques for the systematic analysis of power-corrections in the 
full four parton phase-space are not well established yet. The discussion above clearly 
shows that such techniques are necessary for the analysis of event-shape variables. 
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t^The shape- function approach resums the power-corrections to the distribution which are the 
most singular in the smaU t hmit, i.e. corrections that scale as l/((5"i") for any n. Less singular 
corrections to the distribution, e.g. l/((5"i"^^) can be quite important for the moments. 
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A. Asymptotic expansions for the characteristic functions 



The asymptotic expansions of JF<im>(e), calculated according to eq. ( P78| ) or ( |3.!j| ) 
are the given in eq. ( [AID through ( [OD below. Table | summarizes the numerical 
values of the coefficients. 



<t>lej 



-2 dilogS — - TT' 



In 3 

6 36 8 

4v^+[2 + 61n 3] e-?e=^/2 



(A.l) 



9 



+ 

+ 



28 



— In 2 + 2 dilog3 - -In - +77^ + 7^ 



17 1 



14 



4 



- In ( - ) + - In 2 - - 
1 



12 6 
64 



e2_4e5/2 



11 



105 



In 2 In - + 



2719 
1260, 



(e) = _2 dilog 3 - - TT^ - - In 3 + — 

<t >v ; ^ 6 8 216 

91 127 
12dilog3 + 7r2-- + — ln3 



(A.2) 



+ 
+ 
+ 
+ 

•^<t3>(e) 



16 



-3/2 



11 



TT 



1 In 



4 ve 

8dilog3 + — In 2- - In 
^15 3 



22dilog3 - 13 In 3 



104 , ^ 4 
In 2 + - 

3 3 



16 



,5/2 



12 



In 



1 

e 

1003 152 , ■ 
+ — ln2 



103 



30 



TT 



315 
1 



105 



6^ + 



+ 
+ 
+ 
+ 



- 9 28 83 , 

-2 dilog 3--7r=^ + In 3 

^ 6 135 64 

!^2^1n3-- + 32dilog3 
3 16 9 ^ . 

227 23 

+ — TT^ + 76 In 2 + 46 dilog 3 - — In 3 

12 6 ^ . 

1127 302 



(A.3) 



9 

155 



180 



+ In 2 + 7r^ + 12dilog3 



15 



4 
9 



128 
225 



,5/2 



-- In ( - 



34999 1201 , 

+ In 2 



5040 



105 



1 



1259 649 



(g) = _2 dilog3 - - TT^ + - — In 3 
' ^ 6 4860 480 

85 3373 , , , 2 

^ In 3 + 5 TT^ 

9 80 

[4879 400 5 2 326 „ 
H — ^ m 2 + - TT — m 3 + 30 dilog 3 



+ 
+ 
+ 
+ 



60 dilog 3 
52i 



108 3 

-— - — In 2 - — TT^ + 24 In 3 - 64 dilog 3 
30 15 3 ^ 

5849 372 8 2 _ 32dilog3 

315 35 3 ^ . 



(A.4) 



32 



,5/2 



75 
o 32 
^ +75^' 



24 





iri. 






<~.h > 




t 




-LD 7888 
-|-U. ( ooo 


-LD 071 
-|-U.U ( iO 


-ld ni 1 on 

-hU.Ui iZU 


-ld nnoi Q7 
-hu.uuziy ( 


^1/2 


1 in 


A 

— 4 








^1 
t 




-L8 c^QI 7 


n 7QQQ 


-hU.UOOO ( 


-|-u.uiyuoo 


g3/2 


-L/ Vet; 


-8 8889 


+ 1.7778 


-0.4444 




f2 Infill 




—1 ^ 








^2 
c 




\R fic;7c; 
-rU.uo ( o 




_i_n 7fi777 


— U.ZiOOUO 






-4 


+1.7778 


-0.56889 


+0.426667 


e3 ln(i) 




+2.6667 


-0.3333 










-0.1531 


+0.3399 


+0.32290 


-0.429646 


^7/2 








-0.4444 


+0.426667 



Table 6: Numerical values of the coefficients in the asymptotic expansions of the 
characteristic functions for various moments of i = 1 — T. For example, .^<t>(e) = 
0.7888 -4^6 + 8.5917e + --- 



Substituting only the upper limit y/e in the second integral in eq. (p.9|) we can 
isolate the large-angle soft-gluon contribution to J?-'<im>(e). The results are the fol- 
lowing: 
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Note that (|A.5| ) presents exact results rather than asymptotic expansions at small e. 
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